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Abstract 

A realistic physical model for the so called rate promoting vibration (RPV) at 
enzyme action is constructed. The origin of the RPV is assumed to be an oscillating 
electric field produced by long-lived localized vibrational modes in protein dynamics, 
namely, by the so called discrete breather (DB) in secondary structure. The strength 
of interaction of the RPV with the reaction coordinate is evaluated and its effect 
on the reaction acceleration is assessed within the framework of modern theory for 
thermally activated escape rate at periodic driving. We reveal the phenomenon of 
resonant activation in our model elucidating why the frequency of the RPV in the 
range 100 -j- 200 cm" 1 was chosen by the evolution of enzymes as an optimal one. 
The effect of the RPV on the reaction acceleration is shown to vary from moderate 
one (up to 10 3 10 4 ) in the case of three-site DB to enormous (up to 10 6 10 8 ) in 
the case of five-site DB and thus can significantly contribute into enzyme catalytic 
efficiency. Also the model is shown to be compatible with the known functional 
dependence of enzymatic reaction rates on solvent viscosity. 
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1 Introduction 



Comprehension that specific protein dynamics contributes somehow into enor- 
mous reaction acceleration (up to 10 s 4- 10 15 ) by enzymes has a long history 
(see [1], [2], [3], [4], [5], [6], [7], [8] and refs. therein). It is supposed that an 
enzyme provides more efficient reaction activation compared to corresponding 
nonenzymatic reaction due to a dynamical mechanism and that the latter is 
significant enough to be a rival to the effect of lowering of the potential energy 
barrier (due to transition state stabilization) caused by specific interactions 
of the substrate molecule in an enzyme active site with catalytically active 
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groups. The latter paradigm going back to Pauling is a cornerstone of tra- 
ditional chemical enzymology [9] and is developed in a number of vigorously 
conquering approaches such as transition state stabilization by: i) electrostatic 
field of preorganized dipoles in an enzyme active site [10], [11], ii) low energy 
barrier hydrogen bonds [12], iii) entropic, strain and desolvation effects (lead- 
ing mostly to ground state destabilization) [9], etc. Chemical enzymology is 
a vast field of activity and it will not be touched upon here. Instead we deal 
with physical aspects of enzyme action which are still poorly understood. The 
reason for that seems to be in the lack of adequate theoretical tools to evaluate 
the contribution of a dynamical mechanism in the reaction rate and in meager 
direct experimental evidence for its existence. The latter is altered to some 
extent by recent experiments [13], [14], [15], [16], [17], [18]. In these papers the 
importance of protein dynamics for hydrogen tunneling transfer enzymatic re- 
actions was shown and the notion of the rate promoting vibration (RPV) was 
coined [13], [15] (some authors call it protein promoting vibration). To regret 
hydrogen tunneling transfer is involved only in minor and marginal part of en- 
zymatic reactions and up to now there is no direct experimental evidence for 
the significance of the RPV in covalent bond cleavage between heavy atoms 
that is the most typical rate limiting step at enzyme action. For such en- 
zymes the well known solvent viscosity effect on the reaction rate (see [5], [19] 
and refs. therein) still remains the most vivid testimony for manifesting itself 
of protein dynamics in enzyme catalysis. It is commonly accepted that this 
phenomenon is mediated by protein dynamics but the detailed mechanism is 
obscure. 

An important question is what kind of motion in proteins is responsible for the 
dynamical contribution into enzyme catalytic efficiency? Some authors relate 
protein specific motion producing the RPV with conformational motion of pe- 
culiar amino-acid residues which may belong to enzyme active cite or be at a 
distance from it [17], [18]. In our opinion this scenario encounters difficulties 
at attempts to reconcile it with the experiments [20], [21], [22], [23]. The latter 
unequivocally testify that enzymes retain their catalytic capacity at temper- 
atures below the so-called dynamical "glass" transition where conformational 
motion is considerably suppressed (see, e.g., [24], [21] and refs. therein). Some 
guess on the origin of the RPV can be obtained by considering its frequency. 
The authors of [15] conclude that the estimated dominant peaks in the spec- 
tral densities of the RPV indicate motions on the 150 cm -1 frequency scale. 
There are another assessments of this value in literature: < 300 cm" 1 [6] 
and 200 cm" 1 [25]. In the latter paper the RPV is attributed to the mode of 
Amide- VII of a peptide group vibrations. However this mode requires some 
torsion of the plane of a peptide group and is very energetically consumptive 
and unfavorable. In our paper [26] it was shown that there is another possi- 
ble type of long-lived localized vibrational modes in protein dynamics with 
required frequency, namely, the so called discrete breather (DB) in protein 
secondary structures. The non-linear localized excitations named DBs or else 
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intrinsic localized modes are time periodic spatially localized oscillations with 
significant amplitudes of several units in a chain of weakly coupled non-linear 
oscillators while others are at rest or oscillate with negligible amplitudes. They 
were discovered by Sievers and Takeno in 1988 [27] and proved by MacKay 
and Aubry [28] to be structurally stable. By now they are well understood 
and commonly appreciated as a generic phenomenon in nature (see [29] and 
refs. therein). DBs have become a new and very fruitful paradigm in nonlin- 
ear physics and in particular much work is being carried out at present to 
construct models for DBs in biomolecules. The paper [26] is a development 
along this line. It is shown there that the frequency of a DB in an a-helix can 
be obtained equal to 115 cm -1 in accordance with the results of the experi- 
ments of [30] on far infra-red laser pulse spectroscopy of proteins. The DB is 
actually an oscillation of the planes of some peptide groups in an a-helix or 
/9-sheet around their equilibrium positions with considerable amplitude (up to 
10° -i- 15°) while neighboring peptide groups oscillate with much less ampli- 
tudes and more distant groups oscillate with negligible ones or stay at rest. 
Thus the DB is assumed to store and utilize the energy released at substrate 
binding by an enzyme. It is created by an external cause (binding energy) 
rather than by equilibrium thermal fluctuations though can be fed by the lat- 
ter. As a peptide group is known to have a large dipole moment (~ 3.6 D) 
parallel to its plane [31] the participation of the latter in the DB creates an 
oscillating electric field. Such fields (the so called electrostatic fluctuations) 
are for long supposed to be a key factor for enzyme catalysis [4]. The field 
can interact with the reaction coordinate because during the movement of 
the system along the latter a separation of charges takes place yielding the 
dipole moment of the reaction coordinate. Thus protein dynamics affects the 
potential surface for the reaction making the latter to be non-stationary. The 
shape (e.g., the height) of the reaction energy barrier becomes time depen- 
dent. As will be discussed below the oscillating electric field in accordance 
with the general theory of the effect of periodic driving force [32] ""heats up" 
the system by changing its effective temperature thus giving rise to lowering 
of the activation energy of escape which can be much bigger than the real 
temperature even for comparatively weak fields". 

The aim of the present paper is to construct a realistic physical model for 
the dynamical mechanism contributing into enzyme action and to evaluate its 
possible catalytic efficiency. We attain the latter within the framework of the 
archetypical model for the reaction rate, namely that of the overdamped limit 
of the Fokker-Planck equation (see [33] and refs therein) with non-stationary 
force field. This equation is a standard tool for studying phenomena with time 
dependent potentials mostly within the context of stochastic resonance [34], 
[35], thermally activated escape problem at periodic driving in the classical 
regime [36], [37], [38], [39], [40], [32], [41], [42], [43], [44], [45], [46] and resonant 
activation (see [47] and refs therein). For the estimate of the escape rate we 
apply the results of the papers [37], [39], [46]. 
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The paper is organized as follows. In Sec. 2 we consider the oscillating electric 
field created by a DB in protein secondary structure. In Sec. 3 we evaluate the 
ability of this field to accelerate an enzymatic reaction. In Sec. 4 we argue that 
our model is able to shed new light on the origin of solvent viscosity effect on 
enzymatic reaction rates. In Sec. 5 the results are discussed and the conclusions 
are summarized. 



2 Physical model for the rate promoting vibration 

The possibility for existence of a stable zig-zag shape DB in protein secondary 
structures (a-helix or /5-sheet) was shown in the paper [26]. The DB in our 
case is actually the out of phase oscillations of the planes of some adjacent pep- 
tide groups (due to rotations round their torsional angles) with considerable 
amplitudes up to ~ 10° -j- 15° from their equilibrium values while neighboring 
peptide groups oscillate with much less amplitudes and more distant groups 
oscillate with negligible ones or stay at rest. The dimensionless form of the 
equation of motion for the small deviation Xi(t) (x < 0.3) of the angle of 
the plane for the i-th peptide group from its equilibrium value in a chain 
i = 1, N is 



where p is the coupling constant, U{xj) is a nonlinear local potential and 
the dash denotes a derivative in coordinate. In [26] the following results were 
obtained: 

a-helix: p = -0.00193; U(x) « 0.5a; 2 + 0.017a; 3 + .... 

anti-parallel /5-sheet: p = -0.00115; U(x) « 0.5a; 2 + 0.004a; 3 + .... 

parallel /5-sheet: p = 0.0028; U(x) « 0.5a; 2 - 0.3a; 3 + .... 

Since parallel /5-sheet is practically not represented in natural proteins [48] we 
omit it from further consideration. The linearizeation of (1) yields the disper- 



sion relationship ujk(p) = y 1 + 2 \p\ cos(k) where k is the wave number. The 
DB frequency and its harmonics must not coincide with linear modes spectrum 
[29], i.e., satisfy the requirement u DB > u k (p). There is no analytical solution 
of (1) for the general case of a multi-site DB [29]. However numerical simula- 
tions show that in our case of the DB in protein secondary structure the time 
dependence Xi(t) is qualitatively rather similar to harmonic oscillations [26]. 
The latter can be easily anticipated from (1) taking into account the small- 
ness of the coupling constants p and of the anharmonic corrections in the local 



p(xi-x + x i+1 ) - U'(xi) 



(1) 



dt 2 
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potentials. That is why we restrict ourselves by a harmonic approximation as 
sufficient for further analysis 

Xi (t) = x?f(uj DB t + (f>i); fit) = sin(t) (2) 



where x™ is the amplitude and 0; is the phase. At first we consider the simplest 
case of a three-site DB. In this case we can take into account only the motion of 
the central peptide group in the three-site DB undergoing maximal deviations 
and neglect others because in such DB the amplitude of the central peptide 
group by an order of magnitude exceeds those of the neighbors x m = x™ >> 
x™ tl [26]. Returning to dimensional time we should replace ujdb by dimensional 

frequency f2 (f2 = ujdb\J — ef / ) where V e ff is the effective local potential and 
I is the moment of inertia of the peptide group. In [26] the following results for 

y -^j — were obtained : 114.75 cmT x for a-helix; 84.34 cmT x for anti-parallel (3- 

sheet and 85.49 cm" 1 for parallel /5-sheet, i.e., approximately 2 x 10 13 s _1 for all 
of them. Now let us take into account that the peptide group is known to have 
a large dipole moment (p ~ 3.6 D) parallel to its plane [31]. In an a-helix the 
dipole moment of a peptide group is approximately parallel to its direction. If 
we place the origin of the Cartesian frame in the dipole moment and direct the 
axis X along the dipole moment at equilibrium position of the peptide group 
(and consequently along the axis of the a-helix) then the dipole moment at 
time t has the projections: pf{t) = p cos(x m f (Qt)) , pj(t) = p sin(x m f (Qt)) , 
pf(t) ~ 0) (see Fig.l). Taking into account that x m f(flt) « 1 we have: 
pi(t) ^p e x +p x m f{ttt) e y . 

We are going to examine the influence of the electric field produced by this 
oscillating dipole moment on the enzymatic reaction, i.e., on the covalent 
bond breaking between heavy atoms that is the most typical rate limiting 
step at enzyme action. We assume that the reaction coordinate q is located 
at a distance Ri from the dipole moment in the point Ri = Ri e> and its 
direction constitutes the angle a with the axis Y, i.e., (e q ■ e>) = cos a (see 
Fig.l). The quasistationary electric field (Q << c/Ri where c is the speed of 
light) from the dipole moment in this point is 

3(p t (t)-R t ) 11,-Rfa p _ 2 p x m f(Qt) _ 

m = 7m = 7m ex + ^m — ey (3) 



where e is the dielectric constant for protein interior (e ~ 3.5 [31]). We at- 
tribute the dipole moment d = e q e q to the bond to be broken at the reaction 
where e is the partial charge. The energy of the interaction of the electric field 
with the dipole moment d is 

H int (t) = -(d- E(tj) = const - 2 e P xm cos a q /(rw) . (4) 



5 



The factor before qf(Vtt) in the right hand side is the strength of the interac- 
tion. 

Now we can qualitatively consider the case of a multi-site DB where several 
peptide groups oscillate with non-negligible amplitudes. For instance in a five- 
site DB three of them have significant values compared with the neighbors 
(x™ ~ >> Xi± 2 )- This example is quite representative and we further 
restrict ourselves by it. On the one hand peptide groups in the zig-zag DB un- 
dergo out of phase oscillations, i.e., the deviations of the angular variables have 
opposite signs (sgn(Sxi(t)) = —sgn(6xi±i(t))). On the other hand the struc- 
ture of an a-helix is such that there are three peptide groups in its perpendic- 
ular cros-section, i.e., the radial angle between adjacent groups is 120° > 90° 
(actually there are 3.6 peptide groups in the perpendicular cros-section, so 
that the angle is 100° [31] but this subtliety does not modify our futher 
qualitative conclusions). Two circumstances compensate each other. That is 
why at opposite deviations of the angular variables the projections of linear 
deviations of the dipole moments of adjacent peptide groups on the axis Y 
have the same sign, i.e., are in phase (see Fig.l). If we neglect for simplicity 
the differences both in the distances from the dipoles to the location of the 
reaction coordinate and in angles between the dipoles and the direction of 
the reaction coordinate then we can assert that the oscillating electric field 
produced by three peptide groups of the five-site DB has the amplitude that 
is approximately thrice of that for one peptide group of the three-site DB and 
consequently the strength of the interaction in (4) should be multiplied by 
the factor 3. Taking into account the above mentioned differences makes it of 
course less than 3 and we use the qualitative factor 2 3 as sufficient for our 
further analysis. 

In the next Sec. we consider the effect of the time dependent part of the 
interaction (4) on the enzymatic reaction. 



3 Efficiency of the rate promoting vibration for reaction accelera- 
tion 

We use the Kramers' model for a chemical reaction as the escape of a Brownian 
particle with the mass m from the well of a metastable potential (see for review 
[33] and refs. therein) in the presence of the additional time dependent term 
defined by (4). For the metastable potential we choose without any serious loss 
of generality the conventionally accepted simple form V(q) = aq 2 /2 — bq 3 /3. 
This potential has a minimum at q m in = 0, a maximum at q m ax — a/b and the 
barrier height V(q max ) — V(q m i n ) = a 3 /(6b 2 ). In the overdamped limit (high 
friction case) the Langevin equation of motion along the reaction coordinate 
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q is 



dq aq bq 2 2 e p x m cos a 
^ dt mm m eR? 



f(Qt) + yj2jk B T/m £(t) 



(5) 



where 7 is the friction coefficient, k B is the Boltzman constant, T is the tem- 
perature and £(t) is the Gaussian white noise (< £(£)£(0) >= 5(t) where S(t) 
is the Dirac 5-function) with zero mean (< £(t) >= 0). We introduce the 
dimensionless variables and parameters 

q = qb/a; t = — — ; Cl = — — (6) 



and denote 



_ 2 e p x m cos a k B T b 2 k^T 

aeR* /a] ~ a3 ~ 6 [V( q max ) - V(q min )\ U 



where A is the dimensionless parameter characterizing the oscillating electric 
field strength and D is the dimensionless parameter characterizing the inten- 
sity of thermal fluctuations (i.e., the temperature of the heat bath) relative 
the barrier height. Their ratio is the most important parameter of the model 

A 2 e p x m cos a 

15 = lm B T a/b - (8) 



Omitting the overlines we obtain 

^ = q - q 3 + Af(nt) + J2-fDm/a ^(t'ym/a). (9) 
dt v 



The corresponding overdamped limit of the Fokker-Planck equation for the 
probability distribution function P(q,t) [33] is 



Before proceed further we should define the range of the parameters for our 
model. To be supported by evidence we take the numerals for a very typical 
and one of the most studied enzymatic reaction catalysed by Subtilisin. This 
enzyme belongs to serine proteases and brakes the bond between the atoms 
C and N in a peptide group of a substrate of protein nature. At physiological 
temperatures where enzymes normally work (k B T 4.2 ■ 10~ 14 CGS) the 
corresponding noncatalysed reaction typically has the rate constant k ~ 1 • 
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10 8 s 1 while the enzymatic reaction has the rate constant (that of the rate 
limiting step) k cat pa 5 • 10 1 [49]. Thus the total catalytic effect is ~ 10 12 . 



Usually we do not know the friction coefficient 7. However assuming the single 
local-system relaxation time •ym/a [33] to be typically < 0.5 -7- 1.5 ■ 10~ 14 s we 
conclude that for the noncatalysed reaction we must have T = pa 10~ 22 . 
We can evaluate from the Arrhenius factor of the Kramers' rate 



K 



exp 



■l/(6£>) 



that in order to obtain such a reaction rate constant in the absence of the 
oscillating electric field (A = 0) we should take D pa 6 . 22 1 fal0 ~ 3.3 ■ 10~ 3 . 
We also assume the distance from the minimum to the barrier to be typically 
a/b pa 0.5 A. Taking the typical values ( e pa 0.5 of the electron charge (4.8 • 
10~ 10 CGS); p = 3.6D = 3.6 • 10~ 18 CGS; x rn < 0.3; e pa 1 4 3.5) and 
assuming that the DB is arranged favorably respective the reaction coordinate 
(a pa 0) we obtain that for R > 2A the ratio is A/D < 10 for the case of the 
three-site DB and can attain the values A/D < 20 4- 30 for the five-site DB. 
It should be stressed that even in the latter case A < 0.1 at our value of 
D pa 3.3 ■ 10~ 3 and satisfies the requirement \A\ < 0.25, i.e., we are in the 
so called subthreshold driving regime. The latter provides that the potential 
surface always has a minimum and a maximum, i.e., the oscillating field is 
small enough not to distort the physical picture of the chemical reaction as 
the Brownian particle escape from the metastable state. We obtain from the 
dimensional frequency Q pa 2 x 10 13 s~ 1 that the dimensionless one (see (6)) 
is T K « H < 0.14- 0.3. It should be stressed that we are very far from the 
stochastic resonance regime defined by the requirement Q = nTx [35] and are 
in the so called semi-adiabatic regime [42]. Thus we define the range for all 
three parameters figuring in (10) 

D^3.3-10" 3 ; O<0.U0.3; A/D < 10 



for the three-site DB and 

D^3.3-10" 3 ; ft < 0.1 4 0.3; A/D < 20 4 30 (12) 



for the five-site DB. The range of values for noise intensity D in (12) is in- 
accessible for numerical solution of (10) because of enormous escape times 
(1/r^ > 10 22 ). That is why we resort to the analytic results of the papers 
[37], [39], [46] for the escape rate T averaged over the period of the driving 
force oscillations 2ir/Q. We compile the results of these papers in a simple 
phenomenological formula for the prefactor F of the escape rate in the pres- 
ence of the harmonic driving force with the intensity A and the frequency Q. 
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For our potential (V(q) = q 2 /2 — g 3 /3 in the dimensionless form) and in the 
semi-adiabatic regime Y K « Q < 1 it takes the form 





r 


F = 










A 


exp 


~D 



2n 



2tt 



d(f> 



A 
8Dtt 



sin(20) 



(13) 



where T K is the Kramers' escape rate (11) and g(£l) = s/i |"^ ■ For A/D » 1 
the integral in (13) can be evaluated by the steepest descent method and yields 



2irAg{n) 



exp 



Af A3n 

dv ~ Ton 



(14) 



The dependencies of the prefactor F on A/D and Q are depicted in Fig. 2 and 
Fig. 3 respectively. They are obtained from (13) by numerical integration. For 
the potential V(q) = -q 2 /2 + q A /A the function g(fl) is g(Q) = ^Jjj^- We 
further use (14) for the estimate of the catalytic efficiency for the dynamical 
contribution (that is actually the prefactor) at enzyme action. 



4 Solvent viscosity effect on the rate of enzymatic reaction 



It is well known that enzymatic reaction rates exhibit unusual dependence on 
solvent viscosity rj [19] 

hat oc (rj/voT^ (15) 



where 770 — 1 °P, < (3 < 1 that is generally of order of 0.4 0.6 and a is the 
cosolvent molecular weight. This dependence is the more pronounced the less 
is the cosolvent molecular weight, i.e., the higher its capacity to penetrate in 
protein interior [19]. The empirical law (15) is usualy verified in the range 1 > 
i]/i]o < 50 4- 100. As is commonly accepted the mechanism of this dependence 
should be somehow mediated by protein dynamics but its origin is still not 
clear. 

Below we argue that the present approach can shed some light on the empirical 
relationship (15). We recall that for both a-helix and anti-parallel /9-sheet the 
local potentials for the DB are of hard type [26] for which the frequency de- 
creases with the decrease of the amplitude of oscillations. As a result we come 
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to the following scenario. The increase of the viscosity leads to the increase 
of the probability for a cosolvent molecule to be located near the oscillating 
peptide group of the DB that restricts the amplitude of the oscillations simply 
by creating steric hindrance. The more the number of cosolvent molecules in 
the vicinity of the DB the less the amplitude of the oscillations. As a result 
of decreasing amplitude the frequency of the DB is decreased. Now we phe- 
nomenologically find the law for this decrease to obtain the required scaling 
of the reaction rate constant with the solvent viscosity (15). 

In order to obtain from (14) the required empirical law (15) we should have 



exp 



c 



' Ag(n) 



(JL) 
\VoJ 



ftp) 



(16) 



where C is a constant. We take into account that A oc x m and that frequency 
of the DB is some function (that is not kwown in the analytical form) of 
the amplitude of oscillations = £l(x m ). In pure solvent without cosolvent 
molecules (77 = i] ) we have some amplitude of the DB x™ such that the 
frequency mathes the requirement of the resonant activation Q(x™) ~ 0.2-=-0.3. 
Then we can determine the constant C and write 



exp 


^y 9 (n(x-(rj))) 


Y exp 




1 


Jx m (rj)g(n(x m (rj))) 







(17) 



This relationship implicitly defines the decrease of the amplitude of oscillations 
x m (i]) with the increase of solvent viscosity 77. Let us demonstrate it explicitly. 
We suppose that we are in the the region 77/770 >> 1 where Q(x m (i])) < Q{x™) 
so that g(Q(x m (r])) S j — > 1. If we denote e = x m (r])/x™ then we obtain from 
(17) 



(^) 2 %( fi(<) ) e ,J-M^)(, (QW)) _ 



(18) 



Taking into account that most probably ^° ss 10 -7- 15 (see next Sec.) we 
see that the solution of this equation is e — > 0, i.e., indeed x m (rf) is sharply 
decreased copared with Xq at large rj/rj . The construction of a physical model 
for solvent viscosity effect on the amplitude of the DB x m (i]) and for the 
dependence of the frequency of the DB on the amplitude of oscillations Q(x m ) 
is beyond the scope of the present paper and will be the subject of a separate 
investigation. 
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Concluding this Sec. the following comment should be done. Let us write the 
full expression for the logarithm of the escape rate obtained from (14) and 

(11) 

InT a: 

const _ _L + ±(g(Q) - ^) + ln(yD) - In^nAg^. (19) 

This formula predicts two effects that at first glance can be detected in the 

experiment. First the factor ln(^/~D^J leads to a deviation of the Arrhenius plot 

from the straight line. However the deviation is negligible because this factor 
is fully overshadowed by a much more significant dependence Second the 
factor j}g($T) suggests that not only the prefactor, but also the exponent in the 
escape rate should scale with the viscosity of the solvent due to the dependence 
of fl on it. However this dependence also seems to be hardly observable because 
the driving amplitude is very small A << 1/6. In the next Sec. it will be 
argued that at present the traditional chemical enzymology leaves no reasons 
to expect the driving amplitude to noise intensity ratio to be very large. Most 
probably we have A/D < 15, i.e., at our value D ps 3.3 • 1CT 3 we obtain 
A < 0.06. To regret this value seems to be too small compared with 1/6 ps 0.17 
for the dependence of the exponent on solvent viscosity to be noticed in the 
experiment. No matter how insignificant the effect is, this specific prediction 
remains the most direct way to verify the theory experimentally. 



5 Discussion 

We construct a realistic physical model in which the origin of the RPV at 
enzyme action is the oscillating electric field produced by dipole moments of 
peptide groups participating in long-lived localized vibrational modes (i.e., 
in a DB) of protein dynamics. We evaluate the magnitude of the oscillating 
electric field from peptide groups in two cases. In the first case the field is 
produced by a single peptide group that is the central one of a three-site DB. 
In this case we conclude that peptide group can provide the driving force 
amplitude to noise intensity ratio A/D up to 10. For this purpose the DB 
should be favorably situated relative the reaction coordinate, i.e., not far and 
at favorable angle (see (8)). In the second case the field is produced by three 
peptide groups of a five-site DB. In this case we can attain A/D up to 20^-30. 

The efficiency of the RPV as a dynamical contribution into enzyme catalytic 
efficiency is evaluated by us within the framework of the existing theories for 
the thermally activated escape in the presence of the driving force [39] , [37] . It 
is interesting to note that revealing the role of driving at activated escape in 
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biological systems is considered by the authors of the paper [32] as "a funda- 
mentally important and most challenging open scientific problem" . We compile 
the results of the theory in a simple phenomenological formula (14) that yields 
the required retardation of the exponential growth of the prefactor with the 
increase of A/D (see Fig. 2) in accordance with the results of [37]. Moreover in 
this formula the phenomenon of the so called "resonant activation" (see Fig.3), 
i.e., the existence of an optimal frequency for reaction acceleration manifests 
itself. The notion of resonant activation is traditionally used for phenomena 
with the potential modulation by random fluctuations [47]. However recent 
numerical simulations of the Langevin equation [50] demonstrate that this 
phenomenon takes place also for deterministic driving and our result repro- 
duces this computational finding. It should be stressed once more that the 
maximum of the frequency curves at Q « 0.2 4- 0.3 in Fig.3 is not related to 
the phenomenon of stochastic resonance (defined by the requirement Q = nT K 
[35]) since the value of T K in our case is negligibly small (r^ 10~ 22 ). 

With the help of the formula (14) we conclude that the RPV produced by a 
three-site DB can provide reaction acceleration up to 10 3 4 10 4 (see Fig. 2 at 
A/D 10 and il m 0.1 4 0.3). If we assume that the RPV is produced by 
a five-site DB we can get A/D up to m 20 4 30. In this case we can obtain 
enormous reaction acceleration up to 10 6 4 10 8 (see Fig. 2 at A/D ^ 20 4 30 
and n ~ 0.1 4- 0.3). Thus the present model has potential capacity to explain 
significant part of the total effect of enzyme efficiency. Moreover it gives an 
answer to the question why the frequency range ~ 100 4 200 cm -1 of the 
RPV was chosen by the evolution of enzymes as an optimal one. This range 
(dimensionless ~ 0.1 4- 0.3) turns out to be favorite because it matches the 
requirement of resonant activation where driving is most efficient for reaction 
acceleration (see Fig.3). 

Also the formula (14) is shown to be compatible with experimental data for 
solvent viscosity effect on enzymatic reaction rates at a physically understand- 
able assumption that the amplitude of oscillations in the DB decreases with 
the increase of solvent viscosity by the relationship (??). The latter actually 
means that with the increase of the number of cosolvent molecules in the 
vicinity of the oscillating peptide group of the DB the amplitude of the os- 
cillations is decreased because cosolvent molecules create steric hindrance for 
oscillations. For realization of this relationship a separate model should be 
constructed. The latter seems to be a feasible task and is planned for future 
work. Thus in our opinion the present approach sheds new light on the origin 
of the dependence of enzymatic reaction rates on solvent viscosity. 

What kind of experiments can be interpreted with the help of the present 
model? In its moderate form (reaction acceleration up to 10 3 4 10 4 by a three- 
site DB) it can be directly applied to explaining the residual catalytic activity 
of mutant enzymes. It is known that even at substitution of all catalytically 
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active groups in the enzyme active site by inactive ones by methods of gene 
engineering (when one should expect no lowering of the potential barrier) the 
mutants still exhibit residual catalytic activity accurately 10 3 -j- 10 4 [49], [48]. 
Also the known phenomenon of catalytic antibodies can be treated. These pro- 
teins are actually not enzymes. They have no active site and no catalytically 
active groups at all [48]. Still they also exhibit catalytic activity accurately 
10 3 -J- 10 4 [48]. The explanation of chemical enzymology is that they neverthe- 
less bind the transition state of substrate molecules preferentially thus stabi- 
lizing it [48]. The present model suggests another explanation, namely, being 
proteins they contain secondary structure in which DBs producing the RPV 
can be excited at substrate binding. Another case if the extreme estimates of 
the present model (reaction acceleration up to 10 6 -j- 10 8 by a five-site DB) 
were necessitated and called for by difficulties of traditional chemical enzymol- 
ogy. Such enormous acceleration constitutes significant part of total catalytic 
efficiency at enzyme action. Then the present model would be applicable to 
interpreting much wider range of kinetic data. We conclude that the present 
model for the dynamical contribution into enzyme catalytic efficiency is able 
to describe various situations from partial reaction acceleration up to major 
part of the total effect and besides suggests a route to revealing the origin of 
solvent viscosity effect on rate constants. Whether or not nature makes use of 
such mechanism remains the subject of further investigations. 
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Fig. 1. Schematic representation of a part of an a-helix and a reaction coordinate 
q at a distance Ri from it. Here pi is the dipole moment of the i-th peptide group 
and d is the dipole moment of the reaction coordinate. 
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Fig. 2. The dependence of the reaction acceleration factor on the driving force 
amplitude to noise intensity ratio obtained from (13). The values of the frequencies 
f2 of the driving force from the down line to the upper one respectively are: 0.15; 
0.2; 0.25; 0.3. The value of the noise intensity is D = 3.3 • 10" 3 . 
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Fig. 3. The dependence of the reaction acceleration factor on the frequency of the 
driving force obtained from (13). The values of the driving amplitude of the periodic 
force to noise intensity ratio from the down line to the upper one respectively are: 
5; 10; 15. The value of the noise intensity is D = 3.3 • 1CT 3 . 
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